Localization Transition in Incommensurate non-Hermitian Systems 
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A class of one-dimensional lattice models with incommen- 
surate complex potential V(9) — 2[XrCos{6) + iXiSin{9)] is 
found to exhibit localization transition at |Ar| -|- |Ai| = 1. 
This transition from extended to localized states manifests in 
the behavior of the complex eigenspectum. In the extended 
phase, states with real eigenenergies have finite measure and 
this measure goes to zero in the localized phase. Furthermore, 
all extended states exhibit real spectrum provided |Ar| > |Ai|. 
Another novel feature of the system is the fact that the imag- 
inary part of the spectrum is sensitive to the boundary con- 
ditions only at the onset to localization. 
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Incommensurate systems such as the Harper equation 
provide an important class of models exhibiting both 
extended and localized states in one dimension (ID). In 
this paper, we study the localization transition in sys- 
tems with competing length scales where the underlying 
potential is complex. The systems under investigation 
are described by the following class of lattice tight bind- 
ing models (tbm), 

ll^n+l + V'n-1 + 2[XrCOs{9n) + iXiSin{9n)]Tpn = Eljjn, (1) 

where 0„ = iiian + a. Here, a is a constant phase factor 
and a is an irrational number, which we choose to be the 
the golden mean. This lattice model describes a system 
where the period of the potential is incommensurate with 
the periodicity of the lattice. For Ai = 0, it reduces to the 
Harper equation which exhibits a localization transition 
at |Ar| — 1. Recently, non-Hermitian systems have been 
the subject of various theoretical &-ra] and experimen- 
tal JTlg] studies. In certain ID random systems, where 
all states for the Hermitian problem are localized, the 
addition of a complex vector potential has been shown 
B to result in the delocalization of states, accompa- 
nied by the eigenvalues becoming complex. The sys- 
tem we investigate exhibits both extended and localized 
states and hence facilitates the study of non-Hermiticity 
on both these phases. By investigating how the non- 
Hermiticity alters the localization transition and as well 
as the eigenenergies in the complex plane, we are able 
to study the correlation between the nature of the eigen- 
spectrum and the transport characteristics of the model. 
It is implicit in the literature [BJ that non-Hermiticity cor- 
responds to dissipation and decoherence and as such we 



arguably explore the impact of these effects on the local- 
ization transition. Two interesting limits that we explore 
in detail are the case of pure imaginary potential A^ = 
and the case where the real and the imaginary parts are 
equal {\r — Xi) which can be described as the strong and 
the weak dissipative limits, respectively. Further, we also 
study the case of non-Hermitian lattice models with real 
spectra. This is an interesting problem, since certain 
complex potentials in quantum mechanics are known to 
have a real spectrum provided the potentials exhibit par- 
ity and time reversal(PT) symmetry Q. Here, we seek 
the criterion for a real spectrum in non-Hermitian lattice 
models exhibiting localization transition. 

We study periodic boundary conditions (PBC) anti- 
periodic boundary conditions (APBC) and free boundary 
conditions to investigate the sensitivity to the boundary 
effects. As expected, only extended states are sensitive 
to the boundary conditions and this sensitivity to the 
boundary effects can be used to distinguish extended and 
localized states. We use Ai^^, which is the real part of 
difference in the eigenenergies between PBC and APBC, 
summed over all states, to distinguish extended and lo- 
calized states: The extended states are characterized by 
finite value of AE^ while in the localized phase lS.Er = 0, 
reflecting its insensitivity to changes at the boundaries. 
Our detailed numerical study, based on sensitivity to var- 
ious boundary conditions and wave functions shows that 
the non-Hermitian system exhibits localization transition 
at |Ar| + I Ail = 1. This implies that the Hermitian and 
the non-Hermitian part of the potential carry the same 
weight in determining the transport characteristics of the 
model. Fig. (1) shows AE",. and A_Ei for weak dissipation 
limit Xr = Xi = A, which is discussed in detail below. A 
rather intriguing result is that the imaginary part of the 
/S.E is sensitive to the boundary effects only at the onset 
to localization transition. This result found to be true 
for other parameter values implies that the life time of 
the metastable system depends upon the boundary con- 
ditions only at the transition point. 

In order to explore the relationship between the lo- 
calization character and the behavior of the eigenener- 
gies, we have extensively studied of the eigenspectrum in 
the two-parameter (A,-,Ai) space. Figs. (2-4) show the 
variation in eigenenergies with the parameters for some 
special cases. These studies suggest that it is the local- 
ized phase of the non-Hermitian lattice model (1) that is 



characterized by a complex spectrum, in contrast to pre- 
vious results B where the spectra becomes complex when 
the localized states become delocalized. In the extended 
phase, the spectrum is real provided \Xr\ > \Xi. Further- 
more, for \Xr\ < \Xi\, the number of extended states with 
real eigenenergies have finite measure. Therefore, the 
fraction of states with real spectrum is finite in extended 
phase and vanishes in localized phase. This provides a 
new order parameter for the localization transition, as 
shown in Fig. (5). 

In contrast with the earlier results M where a non- 
Hermitian vector potential was found to delocalize the 
localized states of the random system, the addition of a 
non-Hermitian potential to the Harper model which ex- 
hibits both extended and localized states does not alter 
the localization character of the system. Furthermore, 
our results associate a complex spectrum with the lo- 
calized states, in marked contrast with the earlier result 
where complex spectrum implied delocalization. This is 
one of the central result of our analysis. It implies that 
the previous results relating complex eigenvalues and de- 
localization must be understood as special to the type 
of the system investigated, namely random system with 
non-Hermitian vector potential, and may not describe 
the generic property of non-Hermitian systems exhibit- 
ing localization. 

We now discuss the two limiting cases: the strong and 
the weak dissipation limit. The first case corresponds 
to pure imaginary potential(Ar = 0). It is interesting 
to note that the model with pure imaginary potential 
exhibits duality very similar to that of the pure Hermitian 
problem [nol. Under the Fourier transformation (FT), 
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the tbm (1) with X = Xi reduces to 



2i 



t>m+l + (Pm^l + JTTC0S( 



'm I'vm 



iE 



(2) 



(3) 



Comparing this with the original model, we obtain that 



E*{X) = -{e{j). 



(4) 



This implies that the real and the imaginary part of the 
eigenenergies are related as Er{X) = \Ei{^). Therefore 
the case of the pure imaginary potential has the interest- 
ing property that the localization transition interchanges 
the real and the imaginary part of the spectrum. At the 
onset to localization (A = 1) the model is self-dual, with 
Er = Ei. Fig. (6) shows eigenenergies at some values 
of the parameter in the extended phase. Due to dual- 
ity these figures also show the spectrum in the localized 
phase with the interchange of Ej. and Ei. At the on- 
set of localization, the spectrum with Er = Ei resembles 
that of the Harper equation. It is intriguing that even 



in this strong dissipative limit, extended states with real 
spectrum have a finite measure. In this lattice model 
with competing length scales, it therefore appears that 
extended states are essential for obtaining real eigenen- 
ergies. This is to be compared to earlier formal results 
where FT symmetry was a key for obtaining real spec- 
trum for complex potentials. 

We next discuss the weak dissipation limit. A,. = A^ = 
A, described by the following tbm. 



V'n+l + V'n-1 + 2Ae*''"l/'n = E^^^ 



(5) 



As shown in figure (1), the system exhibits localization 
transition at A — .5. The localization threshold is half 
of that of the Harper equation because both the real and 
the the imaginary part of the potential contributes to- 
wards localization. In the extended phase, the spectrum 
is found to be real and identical to that of A = limit, 
namely E = 2cos(0„). This explains the constant value 
of /S.Er in the extended phase as seen in figure (1). In 
the localized phase , the eigenenergies are complex and 
appear to be described by the following expression. 



Er = 2cosh{"i)cos{9n) 
Ei — 2sinh{'~f)sin{9n). 



(6) 

(7) 



Here 7 is the inverse localization length of the system 
which is found to be equal to the corresponding value 
for the Harper equation 7 = log{2X). In the limit A —> 
00, eigenenergies lie on a circle. Therefore, the localized 
phase is metastable with the lifetime determined by the 
localization length. 

The FT of the model (5) can be analysed further as 
the transformation (2) reduces tridiagonal matrix to the 
following triangular matrix: 



2cos{dn 



+ 2A0m_i = E(j)„ 



(8) 



With PBC, the eigenvalues of this triangular marix are 
the solution of the following algebraic equation: 



(2A)^ = 1[{E - 2cos{en)) 



(9) 



For A < 0.5, as A^ — > 00, we obtain E = 2cos{6n), which 
is found to be identical with numerically obtained spec- 
trum of the model(5). For A > .5, the FT of the model 
allows real solutions: E = 2cos(6'„) -I- 2A. These real en- 
ergies were not found to be the solutions of the model 
(5). The FT of the model also exhibits complex solu- 
tions. It is easy to see that in A — > cxd limit, the algebric 
equation (9) has a solution where E/{2X) lies on the unit 
circle which is also the solution for the model (5). This 
indicates a deep relationship between the spectra of the 
model (5) and its FT - the details have proven so far 
elusive. 

Another aspect of model (5) is that at the onset to 
localization A = 0.5, the FT of the model describes the 



strong coupling limit of the fluctuations of the Harper 
equation once the exponentially decaying part is factored 
out ||ll|. This also describes the Ising model at the on- 
set to long range order for ii^ = [|3[. This limit has 
been shown to be universal using renormalization meth- 
ods jO] as well as more rigorous analytic tools |l2| . This 
result therefore establishes the multi-fractal character of 
the FT of the wave function at the onset of the localiza- 
tion transition. It should be noted that the eigenspec- 
trum remains continuous at the localization transition in 
contrast with the Harper equation which is character- 
ized by singular-continuous spectrum at the transition. 
Therefore, the A^ = A^ limit of the model (1) provides 
a new class of incommensurate systems where the eigen- 
spectrum is gapless and remains continuous except in the 
localized phase. 

In summary, the localization transition of the Harper 
equation remains unaffected by the non-Hermitian per- 
turbation. Weakly dissipative system are characterized 
by real eigenenergies in the extended phase. As the 
strength of the non-Hermitian potential increases, the 
number of extended states with real eigenenergies de- 
crease approaching zero at the onset to localization. In 
the localized phase, states with complex spectrum have 
full measure. The localized phase is metastable, with 
state lifetimes determined by the localization length. The 
question of real eigenvalues is determined by both the 
transport character of the states as well as the amount of 
dissipation. An interesting result is that the extended 
states with real eigenenergies always have a positive 
measure while localized states with real energies have 
zero measure. Therefore, the measure of real eigenen- 
ergies provides a new order parameter for localization- 
delocalization transition. 

In contrast to the earlier results on localization in non- 
Hermitian systems [Q , the system studied here associates 
a complex spectrum with the localized phase. Further, 
the non-Hermitian potential does not delocalize the lo- 
calized states. The localized phase of the Harper equation 
remains localized for any amount of non-Hermitian per- 
turbation. This may be related with the fact that the 
non-Hermiticity appears in the diagonal part of the po- 
tential while the the non-Hermitian vector potential stud- 
ied earlier [Q affects the off-diagonal part of the lattice 
model. 

One future aspect of study is the extension of these 
results to a classical non-integrable (perhaps a kicked) 
system, thus investigating the effect of dissipation on 
non-integrable systems exhibiting localization. Further- 
more, it should be noted that the equation (1) is the 
fermion representation of isotropic XY spin- 1/2 chain in 
a complex magnetic field which is spatially modulating 
| fi3[ . The consequences of the localization transition with 
complex spectrum on the magnetic properties of the sys- 
tem is another interesting open question. Finally, local- 
ization transition in incommensurate tightbinding lattice 



models corresponds to a transition to strange nonchatoic 
attractors (SNA) ||lj] in quasiperiodically driven maps. 
Therefore, the results of this paper may have important 
implications in the study of SNA's in complex maps. 
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FIG. 1. (a) and (b) respectively show AEr and AEi vs 
X = Xr — Xi for a = 377/610. The AEi appears to be related 
to the derivative of AEr- 



FIG. 2. (a) and (b) respectively show the variation in 
eigenvalues as a function of A = Ar = A,; for PBC for 
a = 34/55. The extended phase of this non-Hermitian system 
exhibits real spectrum. Note that in contrast to the Harper 
equation, there is a bending and merging of levels at the tran- 
sition. Furthermore, unlike Harper equation, the spectrum is 
not symmetric about A = 0. 



FIG. 3. (a) and (b) respectively show the variation in 
eigenvalues as a function of A = Ar for fixed Ai = .25 for PBC 
for a = 34/55. In the extended phase all states have real 
eigenenergies provided Ar > A;. 



FIG. 4. Same as figure 2 and 3 for pure imaginary poten- 
tial Ar = 0. 



FIG. 5. shows the fraction of states with real eigenenergies 
as a function of A for A; = .25 case (a) and Ar = (b). Here 
o — 233/377. The steps seen in this plot are due to finite size 
effects 



FIG. 6. Figure shows the spectrum for differrent values 
of A for a — 55/89. For pure imaginary potential. Three 
verical columns respectively correspond to A = .5, .75, 1: The 
three rows show (a) Er vs Ei (b) Er vs n (c) Ei vs n(sorted 
independently of the real part) 
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